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Abstract: We study flavor violation in a supersymmetric SU(5) grand unification scenario 
in a model-independent way employing mass insertions. We examine how the quark and the 
lepton sector observables restrict sfermion mixings. With a low soft scalar mass, a lepton 
flavor violating process provides a stringent constraint on the flavor structure of right- 
handed down-type squarks. In particular, fi ^ turns out to be highly susceptible to the 
1-3 and 2-3 mixings thereof, due to the radiative correction from the top Yukawa coupling 
to the scalar mass terms of 10. With a higher scalar mass around the optimal value, 
in contrast, the quark sector inputs such as i?-meson mixings and hadron electric dipole 
moment, essentially determine the room for sfermion mixing. We also discuss the recent 
deviation observed in Bs mixing phase, projected sensitivity of forthcoming experiments, 
and ways to maintain the power of leptonic restrictions even after incorporating a solution 
to fix the incorrect quark-lepton mass relations. 
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1. Introduction 

The Large Hadron Collider (LHC) has started finally, which we hope will be the first 
machine to produce supersymmetric particles directly. At this stage, experimental input 
that is still playing a major role in probing the soft supersymmetry breaking sector and 
that will keep doing so even in the LHC era, is the flavor changing neutral current (FCNC) 
and CP violating processes. From this data, one can extract information on the potential 
new sources of flavor and CP violations in the soft supersymmetry breaking terms (see 
e.g. [1] and papers that cite it). A model of supersymmetry breaking/mediation, possibly 
in conjunction with a model of flavor, should be compatible with this information. In 
particular, the past two years have seen new measurements of Bg-Bg mixing, both its size 
[2, 3] and its phase [4, 5, 6] (the latter still with low precision), which provide new important 
restrictions on the mixing between the second and the third families of down-type squarks 
[7, 8, 9, 10]. On the other hand, a new experiment is going to explore the lepton flavor 
violation (LFV) decay mode ^ ej, squeezing its branching ratio down to the level of 
10~^^ [11], two orders of magnitude lower than the current upper bound. Therefore, it can 
be regarded as timely to update an analysis on supersymmetric flavor violation. 
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An interesting option in this style of model-independent analysis is to work with a grand 
unified theory (GUT). We take the SU(5) group for example. Since a single irreducible 
representation contains both quarks and leptons, their fiavor structures are related. This 
enables us to use both quark sector and lepton sector processes to look into a single source of 
flavor violation. It is entertaining to see which observable is supplying a tighter constraint. 
The outcome can serve as a hint concerning which sector has a higher prospect for discovery 
of FCNC mediated by sparticles. For the scalar masses and trilinear couplings to obey the 
GUT symmetry, the scale of supersymmetry breaking mediation should be higher than 
the GUT scale. We suppose that this scale is given by the reduced Planck scale 
Mpi/i/Svi" 2 X lO"*^^ GeV, or very close to it, as is the case in a gravity mediation scenario. 

This work is by no means the first attempt in this direction [12, 13, 14, 15, 16, 17, 18]. 
Most notably, there is a recent article that has performed an analysis in a similar framework 
[16]. Three differences are worth mentioning. First, we use the aforementioned — > 67 
decay mode to constrain the 1-3 and the 2-3 mixings, in addition to r ^ 67 and r ^ /X7 
which were considered in Ref. [16]. This seemingly unrelated process becomes relevant, and 
highly restrictive in some cases, thanks to the radiative correction to the 10 representation 
scalar mass matrix from the top Yukawa coupling and the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing [19] . As a matter of fact, this mechanism has long been known and included 
in many of the preceding model studies [12, 20, 21]. Yet, this is the first instance of taking it 
into account in a model independent analysis allowing for general flavor mixing of sfermions, 
as far as we know. Second, the authors of Ref. [16] assume that the quark and the lepton 
mass eigenstates at the GUT scale are aligned to a high degree. This may or may not 
be the case if a solution is incorporated for flxing the wrong quark-lepton mass relations. 
Especially, the flrst and the second families are subject to unlimited misalignment in general 
[12]. We propose a method to overcome this obstacle to some extent. Third, we elucidate 
the importance of the gaugino to scalar mass ratio as a key parameter governing relative 
strengths of the hadronic and the leptonic flavor violations. We show expansions and 
shrinks of the territory ruled by each of the two sectors. In addition to these reflnements, 
we include remarks concerning the latest hint of anomaly in the mixing phase of the Bg- 
meson [22, 23]. 

This paper is organized as follows. In Section 2, we spell out basics of flavor physics in a 
supersymmetric SU(5) GUT model. Section 3 presents the procedure of numerical analysis 
and the experimental inputs. In Section 4, we exhibit the exclusion plot of each mass 
insertion, and discuss how one can interpret the plot conservatively when the lagrangian 
has non-renormalizable terms for accommodating the flrst and the second family fermion 
masses. This section also has a collection of upper bounds on the sfermion mixings, as 
well as deviations in selected CP asymmetries allowed by the other constraints. With a 
summary, we conclude in Section 5. One can flnd notations of the soft supersymmetry 
breaking terms and the mass insertion parameters in the appendix. 
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2. SU(5) GUT and FCNC 



2.1 GUT relation between squark and slepton mixings 

Let us begin by reviewing basic elements of a supersymmetric SU(5) grand unification 
model, that are relevant to flavor physics. The superpotential has the Yukawa couplings 
and the right-handed neutrino mass terms, 

Wgvt 3 -^eatcdeX'^jT.^'TfH' + V2 \%HaTt''Fjk - \%N^^jaH'' + ^M'^NiNj . (2.1) 



Matter fields in 10 and 5 representations are denoted by T and F, respectively, 5 and 5 
Higgses by H and respectively, and a right-handed neutrino by A^. The indices a, . . . , e 
run over components of the fundamental representation of SU(5), and i,j = 1, 2, 3 indicate 
the family. Obviously, \u and M^v are symmetric matrices while \o and \n are not. The 
above Yukawa couplings, by themselves, predict mass unification of down-type quarks and 
charged leptons at the GUT scale: 

nie = rrid, rrifj, = rus, rUr = mf,. (2.2) 

Among these, the third relation is consistent with measurements at low energies, while the 
first two are not. One way to explain this discrepancy is to make corrections to relatively 
smaller masses by including the following non-renormalizable terms [24]: 



(2.3) 



where S is the adjoint Higgs multiplet responsible for breaking SU(5) down to the Standard 
Model (SM) gauge group. These terms will contribute to the Yukawa couplings of the 
effective theory below the GUT scale, expressed in terms of the SM fields as 

WssM = Q'^YuUHu + Q^YnDHd + L^YEEHa + L^Yj^NHu + \n^MnN, (2.4) 

where the fields denoted by uppercase letters are components of the GUT multiplets, 

Ti^{Q,U,E}^, F,^{D,L}i. (2.5) 

The Yukawa couplings appearing in the superpotential of (2.4) are related to those in (2.1) 
and (2.3) by 

Yu = Xu + ^ 0/i + ^fi + \fi^ , (2.6a) 
Yd = \d-C (^/ii - lh2^ , (2.6b) 

Y^ = XD-d{hi + h2), (2.6c) 
5 



Yj!j=\n + C- hN, (2.6d) 
5 
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where the superscripts S and A denote the symmetric and the antisymmetric part of the 
given matrix, respectively. The small number ^ is defined by 

«-5^-10- (2.7) 

where a is the vacuum expectation value (VEV) of S, expressed as in 

(S) = o-diag(2,2,2,-3, -3). (2.8) 

The contribution from the non-renormalizable terms makes the difference, 

YD-Y^ = Ch2, (2.9) 

and this can account for the first and the second family quark and lepton masses. 

For this purpose, Ref. [16] does not make use of the O(^) corrections, but they rely on 
Georgi-Jarlskog mechanism [25]. Their scenario corresponds to a case in our work where 
the quark and the lepton mass eigenbases coincide, i.e. Ul = Ur = 1 in the formalism 
spelled out below. 

Note that the proton lifetime depends on the structure of non-renormalizable operators 
[26], thereby imposing a restriction on the parameters appearing in (2.3). There are corners 
of the parameter space in conflict with proton decay experiments. The present work is not 
specific to a particular pattern of those terms and is valid provided that they are Planck- 
suppressed. 

In order to discuss flavor violation coming from the sfermion sector, one should fix the 
basis of matter supermultiplets. One can choose a basis of Tj and Fi fields such that 

Yu = V;^YuU^, Yd = Yd, Ye = uIYeU*r, Yj, = uIvIYn, (2.10) 

where the hat on a matrix signifies that the given matrix is diagonal with positive elements 
[27], Vq and Vl are unitary matrices in the standard parametrization [28, 29] each with 
three mixing angles and one phase, and [/q, Ul^ and Ur are general unitary matrices. 
Note that Yjj may not be a symmetric matrix, unlike \u- In this basis where Yd is 
diagonal. Ye may not be diagonalized in general due to the difference (2.9), and it should be 
decomposed into the above form using Ul and Ur. These two unitary matrices describe the 
mismatch between the down-type quark and the charged lepton mass eigenstates, arising 
from breakdown of the Yukawa unification Y^ = Yd which is a consequence of SU(5) at 
the renormalizable level. Since Ul and Ufi are crucial in correlating hadronic and leptonic 
processes, we need to examine their structures. We can estimate the size of an off-diagonal 
element of Ye in the unit of the tau Yukawa coupling. 

Ye -Yd -jh^ ^ ^ 

— ^ = 77 77V ^ ~ COS/3 /l2 , (2-11) 

[Ye]33 m,/(7; COS/3) 

where v ~ 170 GeV is the Higgs VEV. Notice the suppression by the factor cos/3 for 
high tan/3. Assuming that each element of /12 is not larger than 0(1), one can obtain 
approximate magnitudes of 1-3 and 2-3 mixings [12], 

[C^L]3a ~ - cos/? [/l2]3a, [f^L]a3 ~ COS/3 \U lK^'^l<\o.Z-, 

J 3a) 



{Vnha ~ - COS/? {h\U, {UrU ~ cos/3 {Ulhlul\ 



(2.12) 
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for a = 1,2. Note that they are suppressed by cos/3. The other entries of Ul and Ur 
can be of 0(1). Finally, we relate the fields to the down-quark and charged lepton mass 
eigenstates as 

Q = q, U = U^u, E = Ule, D = d, L = ull. (2.13) 
This leads us to the superpotential, 

WssM = q'^[V^Yu]uHu + q^[YD]dHd 

1 (2-14) 

+ f[YE]eHd + F[VIYn]NHu + -N^M^N. 

One can notice that Vq is the CKM matrix at the GUT scale. If M^v is diagonal in this 
basis, one also has Vl = f^pMNS' Otherwise, the lepton mixing matrix receives additional 
rotations for diagonalizing Mjv. 

Let us turn to the soft supersymmetry breaking sector. The SU(5) symmetry relates 
the soft supersymmetry breaking terms of squarks and sleptons in a single GUT multiplet. 
The scalar mass terms are given by 



-Aoft 3F^m|F + rtm^T + F^^m|F + rt^m^'r+-- - , (2.15) 



in which the higher dimensional terms involving S are suppressed by 0(0- In terms of 
these soft mass parameters of the GUT multiplets, one can express the soft scalar mass 
matrices of the SM fields as 



2 2,1^2/ 2* 2 2 2/ 2* 2 , ^ ^ 2/ \ 

ttiq = rurp + —t^ mrp , rriu = rrirp — -t^ rurp , = rrirp -\- mrp , (i.lbaj 



2* 2 I 2 ^ 2/ 2 2 ^ A 2/ /o 

nip) = mp + -k,mp, m^^ = rUp - -t^rUp, (2.16b) 

using (2.5). From these expressions and (2.13), one can see that the mass insertion param- 
eters of down-type squarks and sleptons at the GUT scale are linked by 

5iL = UL6l^RUl + 0{i), (2.17a) 
5'RR = UR5tLUi + O{0. (2.17b) 

We can notice two possible sources of deviation from the naive equalities [14], 

^LL = ^RR, ^RR = ^tL- (2-18) 

One is the higher dimensional terms in (2.15), which makes the O(^) corrections, and the 
other is Ul and Ur, the unitary transformations parametrizing the misalignment between 
the down-type quark and the charged lepton mass eigenstates. The former type of cor- 
rections is negligible compared to the typical size of a scanning mass insertion parameter 
appearing later on. On the other hand, these corrections might be comparable to the 
renormalization group (RG) contribution to (5^^. Unless they are tuned in such a way that 
they cancel out the RG-generated Si^rr, they nevertheless do not undermine the importance 
of /i — > 67 constraint. The latter needs more consideration. Obviously, Ul and Ur depend 
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on /i2 through Ye- If /12 is diagonal in the basis where Yd is diagonal, Ul and Ur are unit 
matrices, and (2.18) becomes a fairly good approximation correlating squark and slepton 
flavor mixings. If /12 is not diagonal, the correlation gets loose, but in many cases, LFV 
processes can still give meaningful restrictions on the down-type squark mixings, thanks 
to the suppression of 1-3 and 2-3 mixings shown in (2.12). Examples of this situation will 
be presented in Section 4.2. 

In a similar way, the GUT symmetry links the scalar trilinear coupling terms of squarks 
and sleptons so that their chirality- flipping mass insertions have the relations, 

= Ul 6il uj, + 0(0 x Ao{Hd)/fq, (2.19) 

where Aq is the overall scale of the A-terms and ffij is the average slepton mass. In what 
follows, we do not use this expression since we will ignore the ^-term contributions to 
flavor violating processes. 

2.2 RG running of scalar masses 

RG running from one scale down to a lower scale generates off-diagonal elements of a 
scalar mass matrix. For our purpose, we need to consider two intervals of scale: from M=k 
to MguT) and from Mqut (via M/j) to MsusY- The former is needed to determine the 
boundary condition to give on the soft supersymmetry breaking terms at the GUT scale, 
and the latter is to connect the given boundary condition with low energy observables. 

First, we think of running between M^, and Mqut- Using one-loop approximation, the 
RG-induced off-diagonal elements can be written as [30, 31] 

^s^T ^ - 7A2 [3A*c;A?; + 2XhXl]{3ml + \Ao\') In (2.20a) 
(47rj^ Mgut 

Agm| - ^[4A1,Ad + AjvA^](3m2 + |^o|') In (2.20b) 

where mo is the scalar mass and Aq is the trilinear scalar coupling. Let us focus on the mass 
matrix of T fields, which feeds into the mixings of left-handed squarks and right-handed 
sleptons. From (2.6a), (2.10), (2.16a), and (2.20a), one can obtain the following form of 
RG contribution to the LL squark mixing at the GUT scale, 

(4)" = ~. ivi^y^. 1" ^ + OiQ. ,2.21) 

The 0{(,) correction in the second term is not necessarily smaller than the first term coming 
from the CKM mixing and the large top quark Yukawa coupling. Neither is it very likely, 
however, that they cancel out leading to a value much smaller than the first term. That is, 
the left-handed squark mixing in the above expression, without the 0{^) correction, can 
be regarded as the minimal value of {Sfj)LL that is expected in a supersymmetric SU(5) 
model with the cutoff at M^,. Let us record the CKM matrix dependence of the above 
minimal mass insertions, 

iSf2)LL ~ VtdVts ~ A^ (4)ll ~ V;,Vtb ~ A^, {6i^)LL ~ VZVtb ~ A^, (2.22) 
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where we also express them as powers of A, sine of the Cabibbo angle. 

Using (2.17b), one can get the right-handed slepton mixing from (2.21). Again, we 
drop the 0{(,) term in (2.17b), assuming that it does not conspire with the first term to 
result in a drastic cancellation. If Uji is an identity matrix, {8\j)RR has the same pattern 
as (2.22). Otherwise, one should take the misalignment into account. As (2.12) shows that 
the 1-3 and 2-3 mixings are suppressed, one can rephrase (2.17b) into 

= [UrU [f^i?]33 + 0(cos'/3 Si^), a, 6 = 1, 2, (2.23) 

where \UB\ab-, the upper-left 2x2 submatrix of Ur^ is approximately unitary. [Supposing 
universal scalar masses at M^, one actually has another term of the form [f7i{]a3 {^f'iiLL [f^ijJss 
where (with an abuse of notation what we here call) (533)ll is given by setting z = j = 3 
in (2.21). In what follows we discard this term although it can be larger than what is kept 
in the above equation. Even if it happens to be no n- negligible, it generically enlarges the 
rate of — > 67, only to reinforce the sensitivity of this LFV channel.] Keeping only the 
powers of A, one can schematically rewrite this as 

(4)i^i^~[t^i^]llA' + [^/R]l2A^ {S^^) rr ^ [U rUX" + [U rUX^ . (2.24) 

The mixing between the first and the second families, described by \UR\ab^ is not particularly 
restricted to be small. There can be small, large, or no mixing. One finds that ((^^3)/?/? is 
generically not much smaller than A^, unless the mixing is fine-tuned in such a way that 
the two terms cancel out in either of (2.24). For example, the mixing angle should be tuned 
between —A it A^ in order to have K^Js)/?/?] ^ A^. 

Next, we should turn to the running below Mqut- Before examining an off-diagonal 
entry of a scalar mass matrix, let us recall the running of a diagonal element since a 
mass insertion parameter is normalized by it. Squark and slepton masses at MgusY are 
approximately related to the GUT scale variables by 

m\{MsvsY) « (1 + 6x) ml (2.25a) 
^ (Msusy) ~ ml (2.25b) 

with the definition of gaugino to scalar (squared) mass ratio, 

X = M^^^/ml (2.26) 

The squark mass increases considerably by the gaugino mass contribution. The slepton 
mass actually receives a small correction from the gaugino mass, but it can be ignored for 
later discussions. These facts will be crucial to understanding parameter dependence of a 
constraint. 

Unless tan/? is extremely high, an off-diagonal element of does not run significantly, 
while running of the left-handed squark mass matrix makes the difference [32] , 

^^Hh - [ViY^VQ],, (3mg + l^oP)ln^. (2.27) 
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Using these facts and (2.25a), we can associate squark mass insertions at MgusY to those 
at A/qut as 



[(>ij)RR[MsvsY) ~ l + 6x ' (2.28aj 

(dijjLL(MsusYj ~ r+6x ' (2.28bj 

with the definition 

= ^Amlh/ml (2.29) 

In a parahel way, one can relate slepton mass insertions at a low scale to those at a 
high scale by 

{5\^)RR{MsmY) ~ {5\j)rr{Mgvt:), (2.30a) 
(<5-,)ll(Msusy) ~ {5\,)ll{Mgvt) + kj, (2.30b) 

using (2.25b) and the definition lij = As[mf]ij /ttiq with the radiative correction to the 
off-diagonal slepton mass matrix entries [33] , 

As[m?]., ^ [VmVL]^, {3ml + \Ao\') In (2.31) 

This estimate is based on the assumption that the right-handed neutrinos are degenerate so 
that they are integrated out at a single scale Mr. If they are not degenerate, it is modified 
to involve mixings, phases, and eigenvalues of (see e.g. [13]). Even in this case, it has 
been shown that one can use the above form of expression by replacing Mr with the largest 
eigenvalue of Mtv, if there is a large hierarchy among the right-handed neutrino masses 
[13]. Unlike the quark sector, we do not yet have much information on the neutrino Yukawa 
couplings. They can be of 0(1) in the case of heavy right-handed neutrinos, or extremely 
small if the neutrino masses are of Dirac type. Even if we suppose that seesaw mechanism 
is working, a vast range of right-handed neutrino mass scale is possible, from around the 
GUT scale down to the weak scale. Although the lepton mixing angles have been measured 
to an extent, they cannot be directly related to the mixing matrix Vl due to the additional 
degrees of freedom in M^r, the Majorana right-handed neutrino mass matrix. Moreover, 
the hierarchy of neutrino masses is unknown yet. As the magnitude of lij in one model 
can greatly differ from another, we choose to drop it in the following analysis. Therefore, 
the results shown later are legitimate only for a scenario where right-handed neutrinos are 
light enough for lij to be negligible in (2.30b). (For a study on a case with a large neutrino 
Yukawa coupling and a specific boundary condition on the soft terms, see e.g. [18, 34, 35].) 

Nevertheless, there are circumstances where one can tell consequences of non-negligible 
lij. Here, we assume that Ul is a unit matrix. This assumption will be relaxed in Sec- 
tion 4.2. If neutrino Yukawa couplings are large, they affect not only the running below, 
but also above MquTi of "if- Thus, ((^'j)LL(AfGUT) is decomposed into two pieces, 

{SIj)ll{Mgvt) - {Sij)LL{M,) + ah,, (2.32) 
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where the first term represents possible flavor non-universality at the reduced Planck scale, 
and the second is the RG contribution with 



a = 



ln(M,/MGUT) 



(2.33) 



What we will do in the following sections is to search for a set of viable values of ((5-j) ll(Mgut)j 
imposing experimental constraints. In terms of the variables in (2.32), we can interpret 
this procedure in two different ways: we fix lij and scan over {6lj) ll{M^:) , or the other way 
around. As an example of the first option, suppose that one studies a neutrino mass model 
in which the neutrino Yukawa matrix is given, but there is a room for flavor mixing in the 
soft supersymmetry breaking terms. In this case, one can easily guess the allowed region of 
(Sij) ll(Mgut) = {Stj)*RR{MGUT) from the one shown in Section 4.1 using (2.30b): shift the 
region by —hj- This method is applicable to a model with non-degenerate right-handed 
neutrinos as well. Regarding the second option, one can imagine a situation where the 
only source of F mixing is the neutrino Yukawa matrix, i.e. {6lj)LL{M^) = 0. Under this 
condition, (2.30b) can be rewritten as 



which relies on the degeneracy of right-handed neutrinos. Obviously, the allowed region of 
{Sij)LL{MQ\jT) is given by shrinking the one in Section 4.1 by the factor a/ {I + a). 

In this subsection, we used one-loop estimates to understand the qualitative behaviors 
of squark and slepton mixings, but we numerically solve RG equations for quantitative 
analysis in the subsequent sections. 

3. How to impose constraints on scalar mixings 
3.1 Scheme 

One popular way to constrain sfermion mixings in a model-independent fashion is to scan 
over one mass insertion parameter at a time, while setting the other parameters to zero. 
The practical reason to assume all but one of the parameters to be zero is that it is 
difficult or impossible to take more than one complex mass insertions as free variables 
and plot the allowed volume. Despite its makeshift motive, this strategy works as long 
as the parameter being swept by itself makes the dominant contribution to the process in 
consideration. However, there are cases where presence of another mass insertion amplifies 
the contribution from the scanned parameter, thereby rendering the constraint from a 
process much tighter. 

A well known example is — > Xsj. For instance, a single ((Jfs) RR insertion contributes 
to this decay via the gluino loop shown in Fig. 1 (a). If one takes into account nonzero 
{S^-^Irl insertion as well, the diagram in Fig. 1 (b) with double insertions can make an 
additional contribution [10, 36], whose amplitude is enhanced by tan/3 relative to the single 
insertion graph due to the chirality flip on the gluino propagator. The reason for including 
the double insertion diagram, namely considering nonzero {S^3)rl in addition to the {S2-i)RR 




(2.34) 
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('^23)flfl (Xd \ ho (:d \ 



,X X, 
SR hr ' ' ^ 



(a) single insertion (b) double insertion 

Figure 1: Gluino loop contributions of ((^23)^^: to B ^ Xs^- 



under inspection, is not only that it can significantly increase the B ^ Xsj branching ratio, 
but also that {S^^)jiL = rnb{A — /i tan/?) /m~ is generically present and therefore it should 
not be ignored. An s — > (i equivalent has been used in the study of e' /ex [37]. 

Another example is Bg-Bg mixing. This process is affected by {^2^rr as well. How- 
ever, the Bg-Bs mixing constraint on {d2s)RR greatly depends on the size of (523)ii [36], 
and therefore it matters what value of the LL insertion we choose when we are focusing 
on the RR mixing. Apart from the simple-minded choice of vanishing LL insertion, one 
option is to set {62^)11 = {^23) RR [''']' which may be expected from a left-right symmetry. 
Another well-motivated value of {623)11 is the one generated by RG running from the scale 
where the boundary condition is given down to the sparticle mass scale [8]. This value is 
shown in a rather obscure form in (2.28b) and (2.21). It comes from the CKM mixing of 
quark Yukawa couplings and is expected even with universal soft supersymmetry breaking 
terms at M=k. It should be reasonable to expect at least this amount of LL insertion, even 
if one allows for general non-universal boundary condition, which is the case in this work. 

In the framework of supersymmetric GUT, the story / n 

\"33)RL 

can be extended in a more interesting way. The afore- (gi )* r ((J' 

mentioned parameter {623)rr is related to {623)ll at / % 

the GUT scale, and it can lead to LFV. An obvious ^0 
decay mode is r — > /i7 [13, 14, 15, 16]. It can serve ^ 
as another constraint on (^23)^^?' under the assump- 
tion of SU(5) grand unification. A less obvious mode is 

u — > 67. Due to the GUT symmetry, the CKM mixing ^. o at j- r i j- ■ 

' J & Figure 2: Ncutralmo loop contri- 

leads to an off-diagonal element of the scalar mass ma- bution to /i 67 with triple mass 
trix of the entire 10 members, while they run from the insertions, 
reduced Planck scale down to the GUT scale. With the 

help of {6[3)rr produced in this way, one can complete a diagram for fi ^ with triple 
mass insertions shown in Fig. 2. This diagram receives mr/m^ enhancement relative to 
the usual chargino loop since it is proportional to (533)_rl [12, 20, 38, 39]. Therefore it can 
give a strong restriction on {62-^)11 and thus on (523)rR- 

The above examples illustrate how much a constraint on a given mass insertion pa- 
rameter can be strengthened due to the presence of another insertion. Then, the question 
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would be what the reasonable default value of a mass matrix element is, while a particular 
mass insertion is being scanned. In this work, we take the following scheme for choosing the 
default value of a soft supersymmetry breaking parameter: by default, the 10 soft scalar 
mass matrix elements are set to the RG-induced values from the top Yukawa coupling 
and the CKM mixing, and the off-diagonal components of 5 mass matrix are set to zero; 
we ignore scalar trilinear couplings supposing that a loop graph arising from a nontrivial 
^-term does not accidently cancel the contributions considered later. 

Using this scheme, we carry out a numerical analysis taking the following steps. From 
the Yukawa couplings and gauge couplings at the weak scale, those at the GUT scale 
are computed by solving the one-loop RG equations. In this process, neutrino Yukawa 
couplings are ignored. After reaching the GUT scale, we move to the basis of q, u, d, I, 
and e such that Yd and Yg are diagonal. We assume that Ul and Ur are identity matrices, 
and thus q, d, I, e are identical to their uppercase counterparts in (2.13). In terms of the 
superpotential parameters, this corresponds to the case where /12 in (2.3) is such that it 
reproduces the observed down-type quark and charged lepton masses, and is diagonal in 
the basis where Yd is diagonal. Consequences of relaxing this assumption will be discussed 
in Section. 4.2. In this super-CKM basis of down-type quarks and charged leptons at the 
GUT scale, we set the soft mass matrix of squarks to the form. 



2 2 
= mo 



/ 1 {6' 



( 1 (^f2)LL [5f^)LL 




1 {5i^)RR , (3.1 



and we determine the slepton soft masses using (2.16) neglecting the 0{^) corrections. The 
other scalar masses including those of Higgses are universally put to mg. The trilinear scalar 
couplings are set to zero. With these boundary conditions given at the GUT scale, the 
one-loop RG evolution of the lagrangian parameters is performed down to the weak scale. 
In order to fill out the mass matrices of scalars, charginos, and neutralinos, we determine 
/i from the electroweak symmetry breaking condition, choosing the positive sign. We have 
numerically checked that changing the sign of /i does not make a substantial difference. 
Then, we have all the sparticle mass matrices needed to calculate flavor and CP violation 
quantities. We do not use mass insertion approximation, but employ mass eigenvalues and 
mixing matrices, thereby taking account of multiple insertion graphs automatically. For a 
quark sector amplitude, we keep only gluino loops, and disregard parametrically suppressed 
corrections from neutralino, chargino, and charged Higgs exchanges. 

Regarding patterns of the mass insertion parameters in (3.1), we consider the four 
cases displayed in Table 1. A parameter indicated as 'free' is a variable to be scanned over, 
and the other three are fixed at the respective specified numbers, according to the policy 
outlined above. Those numbers have been obtained by solving the RG equations for the 
soft scalar mass matrices with universal boundary conditions at the reduced Planck scale 
down to the GUT scale in a supersymmetric SU(5) model with minimal field content [12]. 
In this procedure, we have ignored effects of non-renormalizable operators on the running 
of scalar mass matrices. The size of {5fj)LL depends on mo, and tan/3, where -M1/2 

is the unified gaugino mass at Mqut- This dependence is taken into account in a plot for 
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a different set of input parameters, although the change from the value shown in the table 
is insignificant. 



Fig. 


\iSf2) LL 




\iSf3)LL\ 


l(5f3)LL| 


\iS'l3)RR\ 


1 ('^23)^^1 


4 


4.8 X 10-" 


-5 


1.5 X 10-^ 


7.4 X 10-3 





free 


5 


4.8 X 10" 


5 


1.5 X 10-'^ 


7.4 X 10-3 


free 





6 


4.8 X 10" 


-5 


1.5 X 10-3 


free 








7 


4.8 X 10" 


-5 


free 


7.4 X 10-3 









Table 1: Values of mass insertion parameters to be given as boundary conditions at the GUT scale, 
for the case with mo = 220 GeV, M1/2 = 180 GeV, and tan/3 = 5. The phase of a fixed {Sfj)LL is 
equal to a.Tg{—V^*Vtj), as can be expected from (2.21). The first column points to the plot of each 
free variable. 

3.2 Observables 

We summarize observables from the quark sector and how we use them as constraints, in 
Table 2. 

The mass splittings of and Bg mesons have been measured with high precision. The 
error of AMb^ is 0.8% and that of AMb^ is 0.7%. However, their theoretical prediction from 
short-distance physics is not so precise. The main obstacle stems from /b^Bb^ {/'^^Bbs) 
which enters the hadronic matrix element of B^-B^ (Bg-Bg) mixing, parametrizing long- 
distance QCD effects. The present uncertainty in lattice QCD calculation is around 30% 
(see e.g. [43] and references therein). A popular way to avoid this large uncertainty is to 
take the ratio AMbJAMb^ since the error in {I'b^Bbs) / {f's^BB^) is much smaller. Still, 
the SM prediction of the mass difference ratio has an uncertainty of about 40% due to the 



Observable 


Measured value 


Imposed constraint 




0.507 ±0.004 ps"^ [23] 


0.507 ps-i ± 30% 


sin 2/3 


0.681 ± 0.025 [23] 


2 a 


cos2;3 


> -0.4 [40] 




B{B^ Xdl) 


(3.1±0.9+J]|±0.5) X 10-*^ [41] 


[5 X 10-^ 10-^] 


AMb^ 


17.77 ±0.10 ±0.07 ps"^ [3] 


17.77 ps-i ± 30% 




-0.30-0.02 PI 


[-1.20,0.06] 


-0.76l0:i, -2.37+°:33 [23] 


[-1.26, -0.13] U [-3.00, -1.88] 


B{B ^ X,7) 


(352 ± 23 ± 9) X 10-6 [23] 


2 a 


q4>k 


0.39 ±0.17 [23] 


2 a 


kxl 


(2.232 ± 0.007) X 10-3 [29] 


|,SUSY| ^ 1 cxp. 


e'/eK 


(1.66 ±0.26) X 10-3 [29] 




\dn\ 


< 6.3 X 10-26 e cm [42] 





Table 2: Constraints from the quark sector on sfcrmion mixing. An empty third column means 
that the second column is used as is. 
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errors in the CKM matrix elements [43]. As a comprehensive way to embrace the above 
uncertainties, we require that each of computed AM^^ and AMb^ fahs within 30% of its 
central value, fixing J'^^Bb^ and /^^Bb^- In spite of the seemingly loose conditions, we 
will find that these requirements play impressive roles, given higher soft scalar mass. The 
uncertainty decreases with the progress of lattice QCD, and is estimated to be reducible 
down to 8-10% with 6-60 tera flops year of computing power [44]. The improved constraint 
from this smaller error is considered as well. 

Although sin 2/3 does not suffer from uncertainty in the AB = 2 matrix element, 
its SM prediction depends on Vub, which has a sizable error. When we require sin 2/? 
to be within the 2 a range of its experimental value, we allow for a 2 o" variation in 
\Vub\ = (4.31 lb 0.30) X 10"'^ [29] as well. As with the magnitude of mixing, we estimate 
effects of an improved measurement of sin 2(3 at a super B factory, under the assumption 
that it will converge to its SM value. We use 2% as a projected error of \Vub\ and 0.005 as a 
stadard deviation of sin 2/3 [45]. Given these smaller errors, the central values of \Vub\ and 
sin 2/3, if they remain as they are now, become inconsistent with each other, reflecting the 
present tension between them. Under this condition, the future sin 2/3 measurement would 
appear to exclude the SM, and therefore it would be hard to evaluate the influence of its 
improved precision. However, there is a claim that the tension can be reconciled within 
the SM [46]. We do not regard this as a signal of new physics, and assume that \Vub\ will 
decrease so that it becomes compatible with the present sin 2/3. 

As for (pBsi phase of Bg-Bs mixing, we try two distinct ways of imposing the 
constraint: (a) using the latest data from D0 at 90% confidence level (CL); (b) employ- 
ing the 90% CL range recently reported by the Heavy Flavor Averaging Group (HFAG). 
Regarding option (b), we choose the one obtained with constraints from flavor-specific Bg 
lifetime and Bg semileptonic asymmetry. (What is denoted by c/ib^ in this work is 
in the notation of HFAG.) We present both of these cases as they lead to very different 
impressions of the results — the D0 range includes the SM prediction of (f>Bs and hence it 
still works as a bound on the room for new physics, while the HFAG range lies outside 
the SM value, thereby indicating the size of extra contribution required to account for the 
discrepancy [47, 48]. In order to compare the power of (pBs measurements at LHCb with 
that of LFV, we suppose that the future central value of (pB^ is given by the SM, despite 
the current hint of new physics at the level around 2 a. We assume that the error of (ps^ 
win be 0.009 at 10 fb"^ [49]. 

Measurement of the inclusive branching fraction B{B ^ ^dl) had not been available 
until its preliminary result was recently reported from BaBar [41]. The precision is still 
low. Considering the experimental and theoretical uncertainties, we take modest upper 
and lower bounds guesstimated from the exclusive branching fraction B{B — > p/oJ^) = 
(1.18 ± 0.17) X 10"^ [23]. Unlike B Xdi, the branching ratio oi B ^ X^-f has been 
measured with a high precision. We impose a 2 o" constraint on it. 

We use QCD factorization [50] to evaluate S^p, the sine term coefficient in the 
time-dependent CP asymmetry of S ^ (j)K [51, 52]. This approach has a source of 
hadronic uncertainty stemming from regularizing a divergent integral in the annihilation 
contribution. We follow the original prescription in Ref. [50], i.e. we replace f^dy/y by 



- 13 - 



Mode Present bound Future bound 

B{fi 67) 1.2 X 10-^1 [54] 10-^3 [11] 

B{t 67) 1.1 X 10"'^ [55] 10^*^ [56], 2 x 10"^ [45] 

B{t /i7) 4.5 X 10-s [57] 10"*^ [56], 2 x 10"^ [45] 

Table 3: Constraints from radiative LFV decay modes. 

Xa = {1 + ge''^)ln{mB/Ah), with Ah = 500 MeV, < £i < 1, and < V3 < 27r [52]. We 
regard Sff^ as consistent with the data if it is less than 2 a away from the central value 
for any g and ip. 

We also incorporate the CP violation parameters ck and e'/cK in Kl — > vrvr in the 
list. Although we do not explicitly scan over a 1-2 mixing, kaon physics can be influenced 
by double or higher order insertions. Imposing those constraints, we require that new 
physics contribution to each does not exceed the measured value in size. Concerning ex, 
its prediction from squark mixing can be made with an uncertainty much smaller than |e^^| . 
However, we are assuming that there may be an arbitrary 1-2 squark mixing, although we 
do not make a plot for it. This is why we are using a rather conservative bound. 

Finally, we examine dn the neutron electric dipole moment (EDM). Recently it has 

been pointed out that this observable can be greatly influenced if there are both LL and 

RR down-type squark mixings at the same time [53]. In order to evaluate dn, we add 

contributions through the down quark EDM, down quark chromoelectric dipole moment 

(CEDM), and strange quark CEDM. 

We use the constraints from the lepton — — — — ; ; 

Observable Measured value 

sector listed in Table 3. The second column 

shows the present 90% CL upper bound on "^CP ^-^^ -'- ^-"^^ [^3] 

each mode. The third column is the prospec- ^CP "O-^^ ± 0.65 ± 0.18 [58] 

five upper bound from future experiments. '-'cP 

The new limit on // — > 67 is the goal of ^cT"' 0-004 ± 0.037 [23] 

MEG at 90% CL. Also, higher sensitivity A^p'''^ 



to r — > 67 and r — > /Li7 is anticipated from 



a super B factory. In Section 4.1, we choose 

c, . , , Table 4: Monitored observables. Precisions at- 

to use 10 ° as the future limit on is (r —> 67 ; • ; ; • 

tamable at a super B lactory are summarized m 

and B{t /X7), between the two numbers Xable 6 
in each row of the table. If one wants to use 

2 X 10~^ instead, the result can be obtained easily: multiply the upper bound on a given 
mass insertion from r ^ 67 or r ^ //7, by 1/ y/E. 

Imposing the conditions enumerated above, we estimate possible deviations in addi- 
tional observables of interest, shown in Table 4. The first three measure time-dependent 
CP asymmetries in radiative B decays. The definition of S^^^ is given by [59], 

^ ^ T(B,{t) ^ K*j) - T{Bd{t) ^ K*^) 

""'^^ r{Bd{t) ^ K*j) + r{Bd{t) ^ K*7) (3.2) 
=A^P cos{AMbJ) + Sep sm{AMB^t). 
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Note that the time-dependent CP asymmetry in — > K*"f in our convention has the sign 
opposite to that in the above reference. This observable is sensitive to a new CP violating 
phase in the right-handed h ^ s transition, such as coming from ((523)r-R- We define S^p in 
a parallel way by replacing K* with p in the expression. This can serve as a 5 — > d analog 
of S^p^ , affected by {Sf^)RR- One might as well use S^'j^^ to investigate {Sf^)Rji, and 
we record its variation. The rest three are direct CP asymmetries in radiative B decays, 
B Xsj, B Xd'y, and B — > Xs+d'y, whose definitions can be figured out by setting 
t = in the above equation. They are complementary to the preceding observables in 
the sense that they can probe left-handed CP violating new physics such as {S2^)ll and 
{^13)11- We quote the measured value of each observable if available. 



4. Results 



4.1 Viable region of each mass insertion 

As a preparation for reading plots of the GUT scale mass insertions, we sketch the process 
amplitudes in terms of these variables. This will help us understand how a figure changes 
as a parameter is modified. Keeping only factors of interest, the LFV decay amplitudes 
can be roughly put in the form, 

N ^tan/?- (43)ll(Mgut) . , ^ x 

A{^^ ^ e7) oc ^^^^^ x fe)Ki.(MGUT) • (^L)ll(Mgut) ^ 
ttIq nig 

where ms is the typical mass of a slepton, chargino, or neutralino in the loop. The first 
factor in the second line is {Si\.^)rl = mr{A — p,tan/3) /"ffi^ rewritten with (2.25b). We 
used (2.30) to replace the other mass insertions by the GUT scale quantities. As to hadronic 
observables, let us pick up AMb^ as an example; other constraints can be understood in 
a similar fashion. The Bg-Bg transition amplitude depends on ('523)aa('^23)bb/''^5 with 
A,B = L, R. For instance, we can use (2.28) to recast one of these combinations at MsusY 
as 

MsuSY ^ + ^^)^^S 



ml 



[{S^3)ll{Mgvt) + • (J|3)fl R(MGUT) 



where ms is the typical mass of a squark or gluino in the loop. Note that qij from (2.29) 
is nearly independent of niQ. Therefore, as we vary mo and M1/2, the scaling property 
of (4.2) is determined by its denominator. The other two combinations with A = B = L, R 
scale in the same way. 

With these ingredients at hand, we begin to interpret the results. First of all, let us 
look at the result of a two-parameter scan to get a taste of two mass insertions. In each 
of the four plots shown in Figs. 3, there are 1-3 and 2-3 mixings with different chiralities. 
The gaugino mass M1/2 at the GUT scale is chosen in such a way that the gluino mass 
becomes 500 GeV at the weak scale. The scalar mass mo at Mqut is set to two different 
values that can elucidate complementarity of the quark and the lepton sector processes. In 
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[(^BktKMouT) 



-0.004 -0.002 0.002 0.004 

l(<5'|'3)Li]('^^GUT) 



(a) mo = 220 GeV, Afi/2 = 180 GeV, tan/3 = 5 (b) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 5 




(c) mo = 220 GeV, M1/2 = 180 GeV, tan/3 = 5 (d) mo = 600 GeV, = 180 GeV, tan/3 = 5 

Figure 3: Constraints on [(a), (b)] the ((^fa)^^, (^fa)/?^.) and [(c), (d)] the ((<5f3)flfl, ('Jsa)^^) 
planes. For each LFV process, the thick curve is the present upper bound and the thin curve is the 
prospective future bound. A Hght gray (yellow) region is allowed by [(a), (b)] AAf^^ or [(c), (d)] 
AMsjj, given 30% uncertainty in the Ai? = 2 matrix element, and a gray (cyan) region is further 
consistent with [(a), (b)] (/f^^ from D0 or [(c), (d)] sin 2/3. The white lines around the center 
mark a possible improved constraint from [(a), (b)] AM^^ or [(c), (d)] AM^^, with 8% hadronic 
uncertainty. The white curves with short thin lines attached to them display a measurement of [(a), 
(b)] (pBs at LHCb or [(c), (d)] sin 2/3 at a super B factory. Those short lines indicate the excluded 
side. Contributions from RG evolution to [(a), (b)] {Sf.^)^^ and [(c), (d)] {823) ll are indicated by 
the vertical and the horizontal hatched strips, respectively. Their widths do not have any meaning. 

the left column, m-o is taken to be 220 GeV, so that the first and the second family right- 
handed down-type squarks have the same mass as the gluino at the weak scale. (The third 
family is slightly lighter.) This is a benchmark case often encountered in the literature on 
supersymmetric flavor violation. In the right column, we change rriQ to 600 GeV. If one 
fixes the 5 parameters at MguT) this mo maximizes gluino loop contribution to B- meson 
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mixing for the gaugino mass chosen here. We elaborate on this point later. On the plots, 
each mass insertion parameter is treated as a real number. 

In each of Figs. 3 (a) and (b), the two axes are {di^)LL = {^isYrr and {Sij'^)]^ = {523)11 
at the GUT scale. Here we restrict the horizontal axis to a range much narrower than 
the vertical axis since we are especially interested in the effect of RG contribution, but 
otherwise the LL mixing can be arbitrary. One can find that {S23)rr is constrained by 
r fij and AMr^ ■ Interestingly, the two plots show different relative significance of these 
two constraints. In Fig. (a), r — > is stronger than AMb^, i.e. a large portion of the 
region allowed by the latter is excluded by the former. In Fig. (b), the order of importance 
appears to be reversed. Although it is early to draw a conclusion since these plots restrict 
the mass insertions to be real, it is obvious that AMr^ gets tighter while t ^ fi^ becomes 
looser if rriQ is changed from 220 GeV to 600 GeV. As rriQ increases, the LFV constraints 
get relaxed because ms gets bigger in (4.1). In fact, ^ in a numerator grows as well, 
but the growth of m| in the denominator wins. In spite of heavier sparticles, hadronic 
constraints get relatively more stringent. To understand why, pay attention to (4.2). Since 
we have fixed the gluino mass while raising itiq, a squark is heavier than a gluino, and 
we should substitute (2.25a) for m|. Then, the gluino loop contribution to Bg-Bs mixing 
scales like x/{l + 6x)^. This factor increases as we decrease x inversely proportional to 
ttiq, unless x is smaller than 1/12, the maximum point. For x < 1/12, the squarks are so 
heavy that they begin to decouple from low-energy processes as in split supersymmetry. 
Note that x = 0.67 on the left plot and x = 0.09 on the right. This explains the narrower 
AMbs band on the right plot.^ Other quark sector processes are enhanced in a similar 
way, as will be shown later. As can be expected from Fig. 2, /i — > 67 restricts the product 
{^iz) LL{52'i) RR^ resulting in the hyperbolas on the plane. Therefore, the /i — > 67 limit 
on {623)rr varies depending on the size of {6f^)LL- If {^f^jLL = 0, {623)rr is free, but 
the restraint grows severer as \{Sf3)LL\ increases. A special case with RG-induced ((5f3)LL) 
marked by the vertical hatched strip, will be detailed shortly. With increasing ttiq, ^ — > 67 
is doubly suppressed by rn^rng in (4.1b). This expands the area within the hyperbola in 
the plots. The width of the AMb^ band is mainly due to the current uncertainty in the 
Bs-Bg mixing matrix element around 30%. The projected bound with 8% uncertainty is 
depicted by the two white lines around {&2'i)RR — 0- -^^S- (a)) it is not yet as tight as the 
present r fij constraint which will be even tighter in the future. However, in Fig. (b), 
the reduced hadronic uncertainty makes the AMr^ bound more restrictive than t ^ fj,^ 
at a super B factory. Nevertheless, it should be kept in mind that r fij becomes more 
sensitive as tan/3 grows while B^-Bg mixing does not. 

We do the same exercise with a different mixture of mass insertions, {Sf^)RR = {^{sTll 
and {S23)ll = {^23)*rr^ to get Figs. 3 (c) and (d). The vertical range is set around the 
magnitude generated by RG running. Here {6f^)RR is bounded by AMr^, sin 2/3, r — > ej, 
and B — > X^'j. In Fig. (c), the current limits from AMb^ and r — > 67 are comparable to 

similar discussion is given in Ref. [34] in a different context, in a more qualitative way. In their 
scenario, the 2-3 squark mixing arises from large neutrino Yukawa couplings. They state that squark loop 
effects on Bs mixing can be more significant for higher mo. However, they do not mention at what point 
of mo this trend stops and squark loops begin to decouple. 
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each other, which are stronger than B X^^. In Fig. (d), AMb^ with the aid of sin 2/3 
leaves a band which is much narrower than that aUowed by r — > 67. Also, the B — > Xd'y 
bound moves inside the t ^ ej bound. These changes, as well as enhancement of the 
other quark sector processes, can be understood in the same way as the difference between 
Figs, (a) and (b). The inner two white lines around {6f^)jiR = indicate the limit from 
AMb^ with 8% uncertainty in the B^-B^ mixing matrix element, and the outer two white 
lines with short thin lines attached to them arise from sin 2/3 at a super B factory. The 
projected limits from AMb^ and r — > 67 are expected to maintain the present tendency of 
relative strengths. That is, they are comparable to each other in Fig. (c) and the former is 
stronger than the latter in Fig. (d). Again, the hyperbolas come from ^ — > 67. Boundaries 
by the neutron EDM appear in Fig. (d) even though all the mixings are real. Extra 
contribution to dn arises from the combination {Sf^) LLiSf^Yjiji of which the LL insertion 
picks up a non-vanishing phase from the CKM matrix, running from the GUT scale down 
to the weak scale. A specific case with RG-induced (^23)^^1 indicated by the horizontal 
hatched strip, will be discussed later. 

Having grasped a picture of how different constraints act on two mass insertions, let 
us examine cases where one of the insertions originates from RG running from the reduced 
Planck scale to the GUT scale. 

We display in Figs. 4, complex versions of Figs. 3 (a) and (b) with the LL inser- 
tions fixed at the numbers shown in Table 1. They indicate regions of {S2^)rr = {S23)*ll 
(in) consistent with observations. This time, tan/? is varied as well as niQ. Let us walk 
through them starting from Fig. 4 (a). A light gray (yellow) area is allowed by AMb^ but 
not by (pBs from D0. A gray (cyan) area is allowed by both. However, most of it is ruled 
out by the LFV processes, as we have already noticed in Fig. 3 (a). One can guess that this 
should be the case even in the near future, comparing the zone surrounded by the white 
curves and the thin circles with their centers at the origin. In particular, the ^ 67 data 
from the MEG experiment should be able to kill all the parameter space except for the tiny 
disk around the origin. It deserves a remark that — > 67 is playing an important role here. 
Being a 2-3 mixing, (^23)-^/? is normally associated with the r — > /X7 process. For example, 
Ref. [16] discusses interplay between leptonic and hadronic constraints in a similar context, 
but they use only r — > /x transitions to restrict {S2^)rr. This difference arises from the 
strategy of setting the mass insertion parameters. Their default value of a mass insertion is 
zero, while our default is the one which is minimally expected from RG running. Therefore, 
they do not find /z ^ 67 limiting a 2-3 mixing as is obvious from Fig. 3 (a). We believe 
that our choice of mass insertions is more reasonable in a scenario where the soft terms are 
generated around the Planck scale such as gravity mediation. It may be argued that the 
RG-induced 10 scalar mixing is not always guaranteed to be sizeable since the cutoff scale 
can happen to be low close to the GUT scale. This is true. However, a low cutoff would 
threaten the validity of making a connection between the quark and the lepton flavors in 
the first place. Non-renormalizable operators shown in (2.3) and even higher order terms, 
generically, give 0(1) contributions to the quark and the lepton Yukawa couplings, thereby 
erasing any trace of their connection in the flavor space as a single GUT multiplet. It 
should be remembered that an RG-induced LL insertion is not only critical to /x — > 67, but 
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(c) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 5 (d) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 10 



Figure 4: Constraints on the complex plane of {S23)rr, with {Sfj)LL generated from RG running 
between the reduced Planck scale and the GUT scale. For each LFV process, the thick circle is 
the present upper bound and the thin circle is the prospective future bound. A light gray (yellow) 
region is allowed by AMb^ , given 30% uncertainty in the AB — 2 matrix element, and a gray (cyan) 
region is further consistent with (pB^ from D0. A thick black curve shows from HFAG. The 
white curves running from top to bottom mark a possible improved constraint from AMb^ with 
8% hadronic uncertainty. The other white lines running from left to right display a measurement 
of (pB, at LHCb. Thin short lines attached to a curve indicate the excluded side. 

also to Bg-Bs mixing. Indeed, the presence of {S2^)ll is rendering the AMb^ constraint on 
{^23)rR tighter [7, 8, 36]. If it were not for {S2^)ll, the gray region would look like the one 
in Fig. 6 (a), where the contribution to AA-Ib^ from (62-1) ll is not enhanced by {S23)rr. 
Another noticeable point is that the information on (pBs from LHCb can play an important 
role in shaping the allowed region. A particular pleasure with this constraint that it does 
not suffer from the hadronic uncertainty that plagues AMb^ ■ Other observables of interest 
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related to 2-3 mixing are S'^p and dn- The area excluded by each of them is depicted. 
Note that dn depends on the phase of ((523)ll as well as on its size. Although the phase 
of the first term in (2.21) is fixed by the CKM matrix elements, the C'(^) correction is 
unknown and may influence the phase of the entire insertion. Varying the phase of {i2z)LL 
amounts to rotating the dn band on the plot around the origin. Finally, the dotted lines 
are contours of S^p . From them, one can read off its largest possible deviation that is 
consistent with the other experimental inputs. Further information on this CP asymmetry 
is collected in Section 4.4. 

Now that we have recognized the general structure of a plot, we try different values of 
parameters. First, tan/3 is doubled from 5 to 10 in Fig. 4 (b). The AM^^ belt does not 
change very much since its dependence on tan/3 is negligible. Each LFV circle halves and 
becomes tighter. This is evident from (4.1), where each decay amplitude is proportional 
to tan/3. The gluino loop diagrams contributing to each of dn and the B cpK decay 
are also proportional to tan/3, and the allowed region shrinks as tan/3 increases. Next, we 
change the scalar mass parameter. In Fig. 4 (c), we raise ttiq to 600 GeV, a value optimized 
for i?-meson mixings. For the reason already explained, the LFV constraints turn weaker 
while the Bg-Bs mixing belt shrinks on the plot. Other quark sector processes are boosted 
as well. Because of this, the impressions of hadronic and leptonic bounds undergo a sea 
change from Fig. (a) to (c). The current // — > 67 limit gets so much relaxed that it does 
not exclude any region compatible with AMb^ and (pBs- One can also notice that dn has 
become much more powerful. Its limit on the imaginary part of {S23)rr is stronger than any 
other bound on the plot. Indeed, the combination of dn and Bg-Bg mixing leaves nothing 
to do for T fj.'j. In the future case, Bg-Bg mixing looks more restrictive than r — > 
and fi — > 67, particularly thanks to improved precision of (pB^ LHCb. This should be 
contrasted with the situation in Fig. (a) where the LFV constraints, both at present and 
in the future, are stronger than the hadronic ones. Lastly, we consider higher mo and tan/3 
in Fig. (d). Each observable changes according to its tan/3 dependence already mentioned. 
For the first time, the B — > Xgj bound becomes visible. Until now, its branching fraction 
has not been sufficiently disturbed by the new physics contribution given by Fig. 1 (b). One 
reason is that this diagram does not interfere with the SM one since they lead to different 
photon helicities. In Fig. 4 (d), however, higher mo and tan/3 cooperate to enhance the 
supersymmetric amplitude. Still, B — > ^^7 is not very restrictive. Its role should be more 
significant for tan/3 much higher than 10. The range of S^p^ predicted in each case is 
summarized in Table 6. 

A remark is in order regarding the recent reports on (pB^ which reveal a small but 
interesting disparity between the combined fit and the SM prediction [22, 23]. Once we 
take the 90% CL range of (j)Bs from HFAG, instead of that from D0, the gray regions 
change to those surrounded by the thick black curves. Since the HFAG result demands 
new physics contribution to 4>Bsi t^is origin on the plane is positioned outside the thick 
black boundary. On the other hand, the mass insertion can be compatible with the LFV 
data only around the origin. This conflict leads to a restriction in an attempt to understand 
the new fit result of 4>Bs with an RR mixing. The trouble is more serious with lower mo 
and/or higher tan/3. As was explained above, lower niQ enhances the LFV branching ratios 
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while suppressing supersymmetric contributions to Bg-Bg mixing. Also, higher tan/3 gives 
rise to higher LFV rates. Figs, (a) and (b) tell us that these cases are disfavored by LEV in 
combination with (j)Bs- The tension between LFV and (/>s^ is relaxed for higher vtiq used in 
the lower plots. Indeed, one can find a small intersection of the (j)Bs area and the r — > ^7 
disk in each of Figs, (c) and (d). This region could be accessed by measurements of r — > ^7 
and /i ^ 67 in the near future. However, the neutron EDM becomes a new obstacle to the 
zone favored hy (pBa as higher mo reinforces hadronic constraints. The limit from dn grows 
more serious for higher tan/?. The problem can be eased by modifying the size and phase 
of {S2'i)LL at the GUT scale to rotate the dn band as already mentioned. 

At this point we stop considering 2-3 mixing of F, and apply the same procedure to 
the 1~3 sector. We present in Figs. 5, complex versions of Figs. 3 (c) and (d) with the 
LL insertions fixed at the values listed in Table 1. They exhibit constraints on {di^)RR = 
i^i3)*LL- start over with Fig. 5 (a). The light gray (yellow) belt is compatible with 
AMb^, which is further reduced by sin 2/3 into the gray (cyan) region. The resulting area is 
completely consistent with B X^j. The width of this area is comparable to the diameter 
of the circle from r — > 67. The restriction from /_f is so strong that it rules out most of 

the gray zone. The fi ^ disk in this plot is smaller than that in Fig. 4 (a). The reason is 
that the decay amplitude is proportional to {S23)rr ~ here, but to {S[^)rr ~ there. 
In a few years, improved lattice QCD should be able to narrow the AM^^ belt down to 
the one between the two white curves, whose width is again comparable to the diameter of 
the future r — > 67 disk. If this narrowed belt is complemented by measurement of sin 2/3 
at a super B factory, the combined constraint could be comparable to or stronger than the 
future r — > 67 bound. The MEG constraint is so tight that the circle appears to be a single 
dot at the origin. The radius of this circle can be looked up in Table 5. The dotted curves 
are contours of S^p. We present its shift that can be expected obeying other constraints in 
Section 4.4. The rest three plots are for the cases with (b) higher tan/3, (c) higher mo, and 
(d) higher mo and tan/3, respectively. They can be understood in the same way as each 
corresponding figure in Figs. 4 was. Let us stress again that with higher mo, the sensitivity 
of hadronic observables to the GUT scale mass insertions is reinforced while that of LFV is 
weakened. In Figs, (c) and (d), the combination of AMb^ and sin 2/3 essentially determines 
the viable areas. This trend is expected to be maintained by a super B factory. One can 
notice that dn in Figs, (c) and (d), and e' /ex in Fig. (d), begin to be visible due to increased 
mo. These quantities are susceptible to the imaginary parts of {Sf^) RRi^is) RLi^ii) ll and 
{^i3)RR{^33)RL{Si2)LL, respectively, although they are not playing important roles here. 
Figs. 5 (b) and (d) show that the amplitude of -B ^ is enhanced by higher tan/3. One 
can see the reason replacing s by d in Fig. 1 (b). 

For the sake of completeness, we report restrictions on the 10 sector mixings as well, 
which can be represented by {5fj)LL = {^[jTrr- In Figs. 6, we examine the 2-3 mixing. 
Comparing them with Figs. 4, one can notice that Bg-Bs mixing is not as restrictive on 
(•^23)^^ ('^23)^?^?- There, the gluino loop contribution from {S2^)rr to the AB = 2 
transition was enhanced by {623)11 from radiative correction. By contrast, {623) ll here is 
not reinforced by (523)-R-R which is set to zero. Nonetheless, AMr^ and cpBs exclude part 
of the plane in Figs. 6 (c) and (d) where their sensitivities to the GUT scale squark mixing 
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(c) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 5 (d) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 10 

Figure 5: Constraints on the complex plane of (5i3)_r_r, with {Sfj)LL generated from RG running 
between the reduced Planck scale and the GUT scale. For each LFV process, the thick circle is 
the present upper bound and the thin circle is the prospective future bound. A light gray (yellow) 
region is allowed by AMb^, given 30% uncertainty in the AB — 2 matrix element, and a gray 
(cyan) region is further consistent with sin 2/3. The white curves running along the belt mark a 
possible improved constraint from AMb^ with 8% hadronic uncertainty. The other white lines 
running across the belt display a measurement of sin 2/3 at a super B factory. Of the two sides of 
a cos 2/3 curve or a white sin 2/3 curve, the excluded one is indicated by thin short lines. 



are maximized. These constraints should be strengthened in the future. The white curves 
with short thin lines attached to them mark an improved (j)Bg measurement at LHCb. 
They appear in all the four cases. The other white curves, appearing in Figs, (c) and 
(d), represent the projected AMb^ limit. Another outstanding point is that B Xs'j is 
excluding larger area of the (523)^^ plane than {S2s)rr. Recall that the supersymmetric 
diagram arising from an LL mixing is added to the SM piece since they have the same 
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(a) mo = 220 GeV, M1/2 = 180 GeV, tan/3 : 
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(c) mo = 600 GeV, Afi/a = 180 GeV, tan/3 = 5 (d) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 10 

Figure 6: Constraints on the complex plane of {523)ll, with (^^2)^^ ^^'^ (<5i3)ll generated from 
RG running between the reduced Planck scale and the GUT scale. For each LFV process, the thick 
circle is the present upper bound and the thin circle is the prospective future bound. A light gray 
(yellow) region is allowed by AMb^, given 30% uncertainty in the AB = 2 matrix element, and a 
gray (cyan) region is further consistent with (f>Bs from D0. A thick black curve shows (pB^ from 
HFAG. The white curves without short thin lines mark a possible improved constraint from AAIb^ 
with 8% hadronic uncertainty. The white curves with short thin lines attached to them display a 
measurement of (f>B at LHCb. Thin short lines attached to a curve indicate the excluded side. 



chirality structure, while they do not interfere in the RR insertion case in Figs. 4. This 
bound grows more stringent for higher tan/3 as is evident from plots in the right column. 
The LFV constraints are noticeably weaker here than in Figs. 4. Concerning r /X7, this 
is because the decay is dominated by neutralino loop here, but by chargino loop there. 
The chargino loop, if present, generically has higher effectiveness per mass insertion size, 
than the neutralino loop. One can find that /.i — > 67 also occurs. It is caused by a 
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neutralino loop graph proportional to {S23)rr{Ssi)rr with RG-induced {S^i)rr- However, 
it is not strengthened by the factor mr/m^, which accounts for the lower branching ratio 
than in Figs. 4. Despite being moderate, the present and future LFV bounds are still 
disallowing portions of the parameter space. The dotted contours show A^^^^ , the direct 
CP asymmetry in S — > Xg^. The numerical value of its variation is shown in Table 4, 
together with that of another related CP asymmetry, A^^^^^'^^^ . 

Now, we switch to the HFAG fit of the phase of Bg-Bs mixing. In Figs. 4 (a) and 
(b), we cannot find a point which falls within the 90% CL range of (j)Bs: even if we allow 
for an squark mixing. Favored regions appear in Figs, (c) and (d), where hadronic 

processes are enhanced. As those regions involve a large 2-3 mixing of left handed down- 
type squarks, they are likely to give a large modification to B ^ Xsj, in particular for 
high tan/?. In Fig. (d), one can notice that a substantial part of the zone of {^2^)11, needed 
to fit (l)BsJ may conflict with S^^. This conflict also grows more serious with increasing 
tan/?. In an attempt to account for the negative value of (pBs with an LL mixing, one could 
have a bigger hope, given a large mixing, higher mo, and low tan/3. Even if this scenario is 
realized, r ^ /i7 and /u — > 67 will be hard to observe even at a super B factory or MEG. 

Finally, we proceed to the exclusion plots on the complex plane of {Sf^)LL = {^isTrr 
in Figs. 7. Let us compare these with those in Figs. 5. A gray (cyan) region here is larger. 
A significant portion of a gray zone is cut out by i? — > X^'j [60]. The LFV circles are 
significantly bigger. Each of the above facts can be explained in a parallel fashion as we 
did in the previous paragraph. The dotted contours are values of A'^^^ . Its discussion 
will follow in a later part. Note that there are cases where the future LFV data may kill 
part of the area that is compatible with B physics measurements. An example is shown in 
Fig. 7 (b) with lower rriQ and higher tan/3. Yet, constraints mostly come from the hadronic 
sector. 

We finish this subsection with a remark on the sizes of the allowed regions shown in the 
preceding figures. We use mass insertion parameters at the GUT scale as the horizontal 
and vertical axes. Therefore, one should be careful in comparing a plot in this paper with 
one from another work, when the latter is using mass insertions at the weak scale. If the 
weak scale variable is a squark mass insertion, one should convert our plots using (2.28) 
beforehand. 

4.2 Non-renormalizable operators and leptonic constraints 

In the numerical analysis of the previous subsection, we have been employing the naive 
relations (2.18). Now, we should discuss how the results will change if we relax this simpli- 
fication and generalize the correlation of mass insertions to (2.17). One could easily guess 
that the one-to-one correspondence between a hadronic and a leptonic channel should be 
disturbed. Yet, it is not completely broken as will be shown below. Tau decay modes 
still limit 1-3 and 2-3 mixings of squarks of either chirality, albeit to a reduced extent. 
Similarly, — > 67 remains a constraint on the RR insertions. 

Let us think about how a tau decay bound should be modified. We first focus on RR 
insertions, and then on LL. If there is an RR mixing, a tau decay amplitude is dominated 
by the chargino loop which is proportional to for r — > 67 or ((^23)^^ ~^ A*7- 
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(c) mo = 600 GeV, = 180 GeV, tan/3 = 5 (d) mo = 600 GeV, M1/2 = 180 GeV, tan/3 = 10 

Figure 7: Constraints on the complex plane of {Sf^LL, with {5^2) ll and ((523)ii generated from 
RG running between the reduced Planck scale and the GUT scale. For each LFV process, the thick 
circle is the present upper bound and the thin circle is the prospective future bound. A light gray 
(yellow) region is allowed by AM^^, given 30% uncertainty in the AS = 2 matrix element, and 
a gray (cyan) region is further consistent with sin 2/3. The white curves without short thin lines 
attached to them mark a possible improved constraint from AA/^^ with 8% hadronic uncertainty. 
The white curves with short thin lines attached to them display a measurement of sin 2/3 at a super 
B factory. Of the two sides of a cos 2/3 curve or a white sin 2/3 curve, the excluded one is indicated 
by the thin short lines. 



Neglecting the 0{^) term in (2.17a), one has 

= [f^L]a6(4)Rij[t^L]33 + 0(cOs'/?4ii), a,h = 1,2, 



(4.3) 



using (2.12), the smallness of mixing of the third family with the other two. [By the same 
token as for (2.23), there can be another term [C/^Jas {^tz)RR [f^LJsa ~ 1-5 x 10~^/ cos/3 -|- 
0.17 cos/3 where we use the (3,3) component of (2.20b) for its estimation. For small 
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neutrino Yukawa couplings, which we assume in the numerical analysis, this term is neg- 
ligible even compared to the smallest upper bound that can be found in Table 5. For 
large Aat, this should be an uncertainty in relating squark and slepton mixings, apart from 
that stemming from running below Mqut-] The mixing between the first and the second 
families, parametrized by [J/^^Jat, is not limited to be small. For instance, consider the case 
where (^23)/?/? is nonzero while {Sf-^)RR is zero, as in Figs. 4. Here, t ^ fi'-f provides a 
significant constraint on (523)-R-R is a unit matrix. Otherwise, it might happen that 

the association of {§2^)11 with {S2'^)rr is weakened by the factor [C/l]22, or in the worst 
case, is completely broken for [C/l]22 = 0. Although r — > /i7 does not occur in this extreme 
situation, r ^ does since (^23)^?^? gives rise to it through {5[^)ll due to the approximate 
unitarity of [Uijab- This argument, for a general [C/lJcj;,, can be summarized in the form, 

\{S{s)ll\' + |(43)llP « \{4,)rr\' + + 0[cos'p{6U% (4.4) 

which determines B{t — > (e + /i)7). The mass insertions appearing above are all at the 
GUT scale. Note that the current experimental bounds on B{t — > fij) and B(t — > ej) 
differ only by a factor of 2.4. Therefore, once one combines these two, one can always give 
an upper bound on each of {62^) rr and {Si^)rr, almost independent of Ul- The error 
caused by non- vanishing 1-3 or 2-3 mixing in Ul, is diminished below 10% even for tan/? 
as low as 3. If one wants to apply this conservative constraint to the case of Figs. 4, the 
radius of each thick t ^ n'j circle should be enlarged by a factor of 1.9. The thick r — > 67 
circles in Figs. 5 should be expanded by a factor of 1.2. Similarly, the future bounds can 
be modified: multiply each by \/2. Even in this case, tau decays remain severe constraints 
on sfermion mixings. 

The same prescription can be applied to the tau decay bound on an LL mixing. Except 
that the amplitude is dominated by a neutralino loop, we can repeat the above line of 
reasoning with L and R exchanged. In this case, a possible additional term in (2.23) arising 
from {S^^)ll, discussed in Section 2.2, is negligible relative to an upper limit from t ^ fij 
or r ^ 67 shown in Figs. 6 and 7. One can obtain a region permitted by r ^ (e + /u) 7 in 
Figs. 6, multiplying the radius of a thick r fij circle by 1.9. The expansion factor for 
Figs. 7 is 1.2. Again, each future bound should be multiplied by \/2- 

Unlike the tau decay modes, /i — > 67 is more involved, and the following method is 
applicable only to an RR insertion. The dominant contribution comes from the triple 
insertion graph in Fig. 2. Including the diagram with opposite chirality structure, we find 
that the decay rate is proportional to 

d ^ l(4)i?fl(4)LLp + |(^)ll(4)/?rI'- (4.5) 
The /i — > 67 data supplies an upper limit on this quantity. One can use (4.4) to show that 

d > mm{\{6[,)RR\M{6^23)RR\'} ' mf3)RR\' + \{43)rr\% (4-6) 

ignoring the term suppressed by cos^/3. In contrast to r — > (6 + /i) 7, /i — > 67 depends on 
the new pivotal factors, {Si^)rr and {S2s)rr- Ignoring the non-renormalizable operators, 
we had their values equal to those of {Sf-^)*ij^ and ((^23)ll Table 1, respectively. As to how 



- 26 - 



{^13)rR {"23)rR change after the non-renormahzable operators are turned on, there 
are three logical possibilities: (a) each value remains at the same order of magnitude; (b) 
either is very small and the other is not; (c) both are vanishingly small. In Case (a), one 
can use (4.6) in order to translate the upper limit on d to those on {6f-^)jifi and ((^23)^2^' 
nearly independent of Ul- We have seen that both {S[r^)fiji and {62'^)rr are at least of 
the same order as {Sf^)LL from (2.24) — otherwise, they should belong to Case (b) or (c). 
Thus, the C/L-independent upper bound on each of (5i3)_r_r and ((^23)^^1 should be given 
by a — > 67 ring in Figs. 4. That is. Figs. 4 are not modified even with this conservative 
interpretation, while the ^ — > 67 circles in Figs. 5 should be replaced by those in Figs. 4. 
In Case (b), the bound inevitably depends on Ul. As above, consider the scenario where 
{623)rr is non-vanishing while {Sf^)RR vanishes. In addition, suppose that (523)-R-Ri foi' 
example, happens to be highly suppressed. Then, (4.3) and (4.5) lead to 

d^\{6{s)RRW6l,)RR\'\[UL]22\'. (4.7) 

The branching ratio scales like |[f7L]22p- Therefore, a ^ 67 circle in Fig. 4 should be 
enlarged by the factor 1/|[C/l]22|- However, we have learned in Section 2.2 that Case (b) 
is not realized unless the mixing angle in [J/fiJab is fine-tuned. In Case (c), which requires 
a conspiracy of Xu, Xd, hi, /12, /i, and /2 in (2.3), as well as the soft terms, fj, ^ does 
not serve as a constraint. 

In the last part of Section 2.2, we discussed consequences of large neutrino Yukawa 
couplings assuming Ul to be an identity matrix. We considered two cases: one where 
neutrino Yukawa couplings are fixed, and the other where boundary condition at M^, is 
fixed at a universal set of values. Here, let us examine how those results change if we relax 
the condition on Ul- For the first case, we include kj into (4.4) to obtain 

\{Sf3)RR\' + \{SIs)rr\' ^ \iS{3)LL - hsl' + |(43)lL " ^23^, (4.8) 

where {Sfj)RR and {SIj)ll are at Mqut and MsusYi respectively. Unless kj is small enough 
compared to the bound on {SIj)ll = {^ij)*RR presented in the previous subsection, the limit 
on the left hand side is appreciably weakened. Note that a model with lij that large is likely 
to be ruled out by LFV data. The second case is more promising. One can extend (2.34) 
in the style of (4.4), to have 

\{Sf3)RR\' + \{SIs)rr\' ^ X [|(4)llP + |(4)llP]. (4.9) 

Therefore, the upper bounds on (523)-R-R ^-^d {6i^)rr attained from (4.4), are further scaled 
down by a/{l + a). 

Recently, an alternative approach to settling down the uncertainties posed by the non- 
renormalizable operators has been reported [17]. Their work in progress makes use of 
dependence of the proton lifetime on the coefficients of the operators [26] in order to find a 
pattern among them. We would say that our strategy is more generic in the sense that it 
relies only on the condition that the non-renormalizable operators are Planck-suppressed, 
although it may not be as predictive as their anticipated outcome. 



-27- 



Mixing 


Fig. 


Present 




Future 






(a) 


7.6 X 10~2-l.4 X 10" 


-1 


1.3 X 10-2 


\{dis)RR{MGVT)\ 


4 


(b) 
(c) 


3.8 X 10-2-7.1 X 10" 
1.7 X 10-1 


-2 


8.1 X 10-3 
4.0 X 10-2 






(d) 


8.3 X 10-2-1.5 X 10" 


-1 


3.9 X 10-2-4.3 X 10-2 






(a) 


2.7 X 10-2-1.4 X 10" 


-1 


2.5 X 10-3-1.2 X 10-2 


\{Sf3)RR{MGVT)\ 


5 


(b) 
(c) 
(d) 


1.6 X 10-2-7.0 X 10" 
4.7 X 10-2 
5.0 X 10-2 


-2 


1.5 X 10-3-7.3 X 10-3 
1.1 X 10-2 
1.1 X 10-2 






(a) 


0(1) 




OiX) 


\{SI,)ll{Mgvt)\ 


b 


(b) 
(c) 
(d) 


0.6-0(1) 
0.7 
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(b) 
(c) 
(d) 


0.6 
0.1 
0.1 




0.1-0.3 
0.06 
0.06 



Table 5: Upper limit on the size of each mass insertion of down-type squarks at the GUT scale. The 
second and third columns indicate the values of mg, M1/2, and tan/3, used in Figs. 4-7. Regarding 
an KR mixing, if there are two numbers separated by a dash, the left one is for C/^ = 1 and the 
right one is for Uh 7^ 1 obeying (2.12). If the two numbers are the same, it is written only once. 
We do the same for an LL mixing on which the alignment condition is given through \Jr instead 
of Ul- For a general Ur, we drop the n ^ ej constraint as we do not have a systematic way to 
impose it. 

4.3 Summary of bounds 

The restrictions on down-type squark mixings at the GUT scale, graphically shown in Sec- 
tion 4.1, are condensed in a numerical form in Table 5. Each number is the maximum 
distance of a point from the origin on the corresponding figure that satisfies all the con- 
straints considered in the present work. As for 4>Bsj ^se the D0 result, which is marked 
in gray (cyan) in Figs. 4 and 6. We would be left with no solution in many cases if we 
used the HFAG fit (which would be a very interesting outcome on its own [34, 35, 47, 48]). 
In order to estimate the power of (pSs measurement at LHCb, we suppose that its future 
central value will coincide with the SM prediction. We make the same supposition about 
sin 2/9. 

Those upper bounds are subject to change of parameters or scheme of uncertainty 
treatment. A variation may also be caused by choosing 2 x 10-^ instead of 10-^ as the reach 
of T ^ 67 and r — > searches at a super B factory. In particular, the strength of a LFV 
constraint depends on Ul and Ur. We take into account the Ul dependence of a maximal 
RR insertion using the method described in the previous subsection. As for fi — > 67, we 
suppose Case (a) therein, i.e. we do not envisage a fine tuning among contributions to 
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{^13)rR °^ {^23)rR- ^ LFV restriction is important, relaxing the assumption of f/^ = 1 
increases the upper hmit of the given insertion. Concerning the hmit on an LL insertion, 
we follow the same procedure to evaluate the dependence of a tau channel onUpt, while we 
keep — > 67 only for Ur = 1. Even if [//j is unity, however, it turns out that the leptonic 
data does not cause a big additional reduction in the bounds set by the hadronic inputs, 
under the conditions considered in this work. A lepton sector constraint should be looser 
if we allow for a different Ur. Therefore, the quoted numbers are not greatly influenced by 
a change of C/r. 

4.4 Possible alterations in observables 

With the region of each mass insertion obtained in Section 4.1, we estimate a possible 
difference of an affected observable from its SM value. The result is summarized in Table 6. 
Four of them have been already displayed as contours on each figure indicated in the table. 
Note that what has been shown as contours is the value of the observable, not the deviation 
from the SM prediction. We use the same set of constraints as in Section 4.3. 

Under the present conditions, there are still CP asymmetries that might potentially 
have a discrepancy bigger than the precision attainable at a super B factory. They are 
Sqp , Sq'p, A^^f^^^ , and A''^^^'^'^'^"' . They show larger possible alterations for higher mo, 
while A^f^^ doesn't follow this tendency. Being hadronic observables, their sensitivity to 
the GUT scale flavor violation is amplified for higher mo, as was explained in Section 4.1, 
although they are more severely restricted by other quark sector processes for the same 
reason. As we did for Table 5, we take account of uncertainties due to a misalignment 
between quarks and leptons of the lighter two families. In this case, we obtain the values 
after the dash signs, which can be larger than the estimates for perfect alignment. 

We repeat the same task with the prospective future inputs. One may expect the pre- 
sented deviations, provided that no constraint is seriously violated in a future experiment. 
With lower mo, S^p^ and S^p will not show a signature detectable at a super B factory, 
even if quark-lepton misalignment is allowed, while A^^^^ and A^^^'^^'^^"' might reveal a 
hint. With higher mo, search for a supersymmetric effect in S^p'^ becomes feasible as well. 

In the case with the {6f^)LL mixing, its effect on A^^^^'^'^^"' is negligible so that the 
variation is at most about 0.5%, because the channel B — > ^5+^7 is dominated by S — > 
Xsj. We include S^^f^^ ^ in the table as well for reference. 

Among the CP asymmetries mentioned above, S^j^ and S(jp are sensitive to RR 
mixings of squarks, and thus are closely related to LFV. Recall that RR mixings give rise 
to much higher LFV rates than LL, as we have seen in Section 4.1. This motives us to look 
into allowed ranges of those two CP asymmetries as functions of LFV branching ratios. 

First, we show the correlation between S(jp/ and B{t ^ ^7) in Figs. 8, each of which 
results from the same set of mass insertions as the corresponding plot in Figs. 4. Every 
point on the figures satisfies the current AM^^ and B{B ^ ^sl) constraints. The upper 
limit on r ^ /i7 from /i — > 67 has been deduced from the contours in Figs. 4. In Figs. 8 (a) 
and (b), what restricts S^p^ at present is r — > /i7, and in the future ^ — > 67 at MEG 
should take over. In Figs, (c) and (d), d^, in addition to r — > /i7, is playing an important 



- 29 - 



Deviation Future 



Mixing 


Fig. 


Present 


Future 


precision 




(a) 


0.04-0.07 


0.007 




\AS^P\ 


(b) 
(c) 


0.04-0.07 
0.18 


0.007 
0.04 


0.02 




(d) 


0.16-0.26 


0.07-0.08 






(a) 


0.06-0.30 


0.006-0.03 




\AS^p\ 

{Sf3)RR 


5 

(c) 


0.06-0.28 
0.21 


0.006-0.03 
0.05 


0.10 




(d) 


0.39 


0.09 






(a) 


0.06-0.28 


0.006-0.03 




{^vi)RR 


5 

(c) 


0.06-0.28 
0.17 


0.006-0.03 
0.03 






(d) 


0.32 


0.05 






(a) 


1.3 


1.3 




\AA^cp^\{%) 


(c) 


1.9-2.3 
3.3 


1.0-1.4 
1.7 


0.4 




(d) 


5.2 


2.8 






(a) 


1.3 


1.3 




\AA}.p^^ (%) 


(b) 
6 ^ ^ 

(c) 


1.8-2.2 
3.2 


0.9-1.3 
1.6 


0.6 




(d) 


5.1 


2.7 






(a) 


16 
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lA^^T,"^^! (%) 
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15 
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Table 6: Maximal departure of each observable from its SM value given the present and the future 
constraints. The second and third columns indicate the plot on which we calculate the observable. 
Of the two deviations separated by a dash in a cell, the left one is for [7^ = 1 and the right one 
is for Ul 1 obeying (2.12), for the first three CP asymmetries. Those two types of deviations 
should be regarded as the same if only one is written. For the rest, the alignment condition is given 
through Ur instead of Ul- For a general Ur, we drop the n ej constraint as we do not have a 
systematic way to impose it. 

role, and the future expectation of S^p'' is determined by AMb^ and (pBs- One can find 
the numerical range of S^p'^ allowed in each of the four figures in Table 6. Note that if 
two numbers are separated by a dash in the table, one should take the left hand side since 
the plots are for Ul = 1- One could translate these plots to a case where Ul is not fixed 
at unity, following the prescription presented in Section 4.2: regard the horizontal axis as 
B{t — > (e + /i) 7)) instead of B{t — > ^'y), and shift the upper bounds on t ^ rightward 
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Figure 8: Correlation between S^p^ and i3(r — > iij) obtained by varying {S23)rr, with {6fj)LL 
generated from RG running between the reduced Planck scale and the GUT scale. A light gray 
(yellow) point is disfavored by neutron EDM, while a gray (orange) point is not, and a black (blue) 
point satisfies the future AMb^ and (pB^ constraints. The dashed horizontal line marks the 2 a 
range of , and its SM value is the solid horizontal line. The present and the future limits on 
T — > /i7 and /i — s- 67 are indicated by the vertical lines. 



in accordance to this change, while keeping the positions of the vertical lines for /x — > 67. 

Second, let us move to the correlation between Sq'p and B{t — > 67), displayed in 
Figs. 9, which correspond to the parameter space considered in Figs. 5. We discard any 
point that is incompatible with the present data of AMb^, sin 2/?, or cos 2/3. The upper 
limit on r — > 67 from /i — > 67 has been inferred as we did in the preceding paragraph. 
For lower mo shown in Figs. 9 (a) and (b), ^ — > 67 provides the limits on S^p both 
currently and in the future. The MEG bound is not visible on the plane since it restricts 
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Figure 9: Correlation between S^p and B{t 67) obtained by varying {Si^)iiR, with {Sfj)LL 
generated from RG running between the reduced Planck scale and the GUT scale. A light gray 
(yellow) point is consistent with all the current constraints, and a black (blue) point satisfies the 
future AMb^ and sin 2/3 constraints. The solid horizontal line marks the SM value of S(7p. The 
present and the future limits on r ^ fi"f and 11 ^ ej are indicated by the vertical lines. In Figs, (a) 
and (b), the MEG line is outside the left border of each plot. 



B(t 67) < 7 X 10^^-*^. For higher mo in Figs, (c) and (d), possible range of S^p is 
determined by the other hadronic observables, with little help from the lepton sector. The 
way to convert these plots to those for Ul ^ 1 is almost the same as above: relabel the 
horizontal axis as B{t (e + /x) 7)) instead of B{t — > 67), and change the upper bounds 
on r — > 67 to those on r — > (e + /x) 7. A difference from the above case is that one should 
also multiply the /i — > 67 limit on B{t 67) by 24 ~ A^^. 

The latest interest in the phase of Bg-Bg mixing leads us to examine its modification 
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0.02-0.03 
0.57 


0.004-0.008 
0.33 




(d) 


0.32 


0.12 



Table 7: Maximal departure of (^b^ from its SM value under the present and the future constraints 
except for those on itself. The second and third columns indicate the relevant plot. Of the two 
deviations separated by a dash in a cell, the left one is for J7l = 1 and the right one is for [/^ 7^ 1 
obeying (2.12), for the RR mixing. Those two types of deviations should be regarded as the same if 
only one is written. In the case with {S23)ll, the alignment condition is given through Uf( instead 
of Ul- For a general Ur, we drop the fi constraint as we do not have a systematic way to 

impose it. 

that can be caused by new physics. We lift the constraint on (pBs while keeping the others 
used in Section 4.3, and record its variation allowed by the other bounds in Table 7. The 
difference between the announced central value and the SM prediction is about 0.7. From 
the table it appears that cases with lower mo and/or large RR mixing (but small LL 
mixing) are disfavored by <j)Bs- the case of RR insertion with higher mo, the primary 
barrier is the neutron EDM as is evident from Figs. 4 (c) and (d). Let us remind the reader 
that this situation can be ameliorated by multiplying {623)11 by an 0{1) complex factor 
at Mqut- With the LL insertion and higher mo, on the other hand. Figs. 6 (c) and (d) 
show that B — > and 5^p exclude a major part of the region preferred by the HFAG 
fit, although there are still remaining parts that are responsible for the large difference in 
(pBs recorded in the table. 

In Figs. 10, we investigate how LFV constrains (pBs by means of correlation plots, 
focusing on the RR insertion case considered in Figs. 4. At first, let us consider only the 
leptonic constraints. In this case, LFV and the latest (pBs fit are better reconciled for 
higher mo depicted in Figs, (c) and (d). Obviously, lower tan/3 is preferable since a LFV 
limit gets tighter for higher tan/3. However, if one takes the neutron EDM bound seriously, 
the light gray (yellow) points are discarded while the gray (orange) points remain, and 
therefore it becomes harder to account for (pB^ with an RR mixing. Remember that one 
can apply this result to a popular benchmark scenario in which the soft terms at M^, are 
flavor-blind and all the right-handed squark mixings are supposed to originate from large 
neutrino Yukawa couplings, as we discussed in the last part of Section 2.2. The recipe is to 
multiply each LFV branching fraction by (1 -|- a)^/a^ with a in (2.33). This factor arises 
from the additional running of slepton masses from Mqut down to M/j, and strengthens 
LFV as the result. 
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One can be more optimistic in viewing the same correlation plots. For example, the 
neutron EDM constraint may be weakened if there is also a non- vanishing complex LL mass 
insertion at M*, or one might simply choose to ignore the constraint due to its hadronic 
uncertainties. Then, it might be that the present status of is hinting at a LFV process 
occurring at a rate that can be explored in the near future. Notice that this scenario works 
best when the value of x defined in (2.26), is around 1/12, as we discussed in Section 4.1. 
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5. Conclusions 



We imposed hadronic and leptonic constraints on sfermion mixing in a class of supersym- 
metric models with SU(5) grand unification. We did not particularly assume that the 

sfermion mass matrices have a universal form at any scale, but rather that any off-diagonal 
entry may be nonzero, which is generically the case in gravity mediated supersymmetry 
breaking. Those off-diagonal elements are encoded in the dimensionless mass insertion 
parameters in terms of which we express experimental bounds on flavor non-universality 
at the GUT scale. While fixing the gluino mass to 500 GeV at the weak scale, we tried two 
different boundary conditions on the diagonal components of the soft scalar mass matrix 
at Mqut: lower mo = 220 GeV and higher uiq = 600 GeV. We varied tan/3 from 5 to 10 
as well. For lower niQ, we have found that the upper limit on an RR mixing is essentially 
determined by a LFV decay mode both at present and in the near future. This is true 
even when one introduces non-renormalizable terms to accommodate the lighter down-type 
quark and charged lepton masses. In particular, the apparently unrelated mode — > 67 
turns out to be remarkably sensitive to a mixing involving the third family. This sensi- 
tivity will be much higher with the progress of the MEG experiment. For higher tuq, the 
situation turns the other way around so that the hadronic constraints, such as _B-mcson 
mixing and neutron EDM, dominate. Also in the near future, measurements at the LHCb 
and a super B factory, with the aid of improved lattice QCD, should be able to probe 
an RR mixing, with a sensitivity higher than that of a LFV experiment. Concerning the 
LL mixings, they are mostly restricted by hadronic data from B physics, although LFV 
supplies additional information if niQ is low and tan/? is high. These findings unveil a nice 
complementarity of the quark and the lepton sector processes showing their strengths and 
weaknesses, depending on the gaugino to scalar mass ratio. We included discussions on the 
consequences of the discrepancy recently observed in the S^-meson mixing phase. 
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A. Notations 

The scalar mass terms in the soft supersymmetry breaking sector of the minimal super- 
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symmetric standard model are given by 

-£soft D Q'^ml Q + U^mfj U* + E^m| E* + o'^ml D* + L^ml L, (A.l) 

where the uppercase letters denote the scalar components of the SM superfields embedded 
in T and F as in (2.5). Consider a basis where the down- type quark and the charged 
lepton Yukawa matrices are diagonalized by superfield rotations. The scalars in this basis, 
denoted by lowercase letters, are related the above fields by (2.13). Therefore, their mass 
matrices are connected to those above by the basis change. 



ml = mQ, ml = UQmfjU^Q, ml = URm\u\f, w?^ = m\), = ULm\ul. (A.2) 

Suppose that the squark and slepton mass terms are given by, 

-C D dsj + e\. [m^ABh ^Bj , (A.3) 

in the basis where the down-type quark and the charged lepton mass matrices are diagonal. 
The sfermion mass matrices include contributions from the Yukawa couplings, the fi term, 
the D terms, the soft scalar mass terms, and the A terms. In terms of the mass matrices, 
mass insertion parameters are defined by [30] 

{6fj)AB = [m|^^]ij/m|, {61^)ab = [m-Ashj/mf, (A ^) ^ (B, j), 

where A, B = L, R denote the chiralities, i,j = 1, 2, 3 are the family indices, and m~ and 
m2 are the average sfermion masses [1]. In this work, we heavily rely on the mass insertion 
notation defined above to discuss the flavor structure of squarks and sleptons. Yet, we do 
not use mass insertion approximation to compute physical amplitudes, but work with mass 
eigenstates and mixing matrices. 

Normally, as its name implies, a mass insertion is a quantity that should be defined at 
the scale of the particle mass. Therefore, a squark or a slepton mass insertion is considered 
at the sparticle mass scale or at the weak scale. This is the case in the previous paragraph. 
In this work, we borrow this notation to deal with the scalar mass matrices at the GUT 
scale: a GUT scale mass insertion is an off-diagonal entry of a soft scalar mass matrix 
divided by the averaged diagonal element, in the basis where the Yukawa matrix is diagonal. 
Following this definition, we have 

{5\j)LL = [mf]ij/mf, {SIj)rr = [mljij/mf, 

for i ^ j. 
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